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Abstract. We report the detection of a thick disk in the edge-on, low surface brightness (LSB), late- type spiral 
ESO 342-G017, based on ultra-deep images in the V and R bands obtained with the VLT Test Camera during 
Science Verification on UTl. All steps in the reduction procedure are fully described, which, together with an 
extensive analysis of systematic and statistic uncertainties, has resulted in surface brightness photometry that is 
reliable for the detection of faint extended structure to a level of V — 27.5 and R = 28.5 mag/sq arcsec. The 
faint light apparent in these deep images is well-modeled by a thick exponential disk with an intrinsic scale height 
about 2.5 times that of the thin disk, and a comparable or somewhat larger scale length. Deprojection including 
the effects of inclination and convolution with the PSF allow us to estimate that the thick disk contributes 20-40% 
of the total (old) stellar disk luminosity of ESO 342-G017. To our knowledge, this is the first detection of a thick 
disk in an LSB galaxy, which are generally thought to be rather unevolved compared to higher surface brightness 
galaxies. 
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^ 1. Introduction 



Outside our own Galaxy, most of what we know about 
the structure, evolution and dynamics of stellar popula- 
tions, and their connection to dark matter, is deduced 
from high surface brightness features: bars, bulges, and 
thin disks. Fainter surface brightness components such 
as stellar halos, thick disks, and globular clusters probe 
galactic potentials differently, in both time and space ow- 
ing to their larger age and extent. The formation mech- 
anisms of these faint tracers are still a matter of some 
controversy; suggestions range from early protogalactic 
collapse, secular processes such as heating from molec- 
ular clouds, black holes and spiral structure, through to 
later stochastic processes such as accretion (see recent re- 



kinematical, morphological and chemical characteristics, 
but too few systems have been sufficiently well studied to 
constrain the models. Due to the difficulty in detecting 
low surface brightness features reliably in external galax- 
ies, the important complementary information they con- 
tain has only begun to be tapped. 

In the Milky Way, faint disk and halo compo- 
nents can be separated on the basis of their kinemat- 
ics and morphology, and - to a certain extent - metal- 
licity, because individual stars can be resolved. The 
Galactic stellar halo of field stars and the globular clus- 
ter systems have volume densities that decrease with 
galactocentric radius r roughly as p(r) 



-3.5 



(Harris fc Racine 1979: Saha 1985: Zinn 1985), sim- 



views by puser 200q ; [Bland-Hawthorn fc Freeman 200q ; jj^r to results for halo populations in large spirals like 



and references therein). These scenarios predict different 
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*Based on observations collected at the European Southern 
Observatory, Paranal, Chile (VLT-UTl Science Verification 
Program) 



M31 ([Racine 1991|; jRcitzcl ct al. 199^ ) and NGC 4565 
(Fleming et al. 1995). Giant ellipticals and supcrluminous 
CD galaxies, on the other hand, which arc thought to 
be the product of many mergers, have halo luminosi- 
ties and globular cluster systems that fall less steeply, 
roughly as p{r) (x r^^-^ (Harris 198(:; Bridges et al. 1991 
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Harris et al. 1995; Graham et al. 1996). The total mass 



r beyond a few 
a review). 



the bulk of which is believed to be contained in dark mat- 
ter halos, is inferred from kinematical studies to have vol 
ume densities that decline as p{r 
disk scale lengths (see packett 1996| for 

Our Galaxy also has a faint thick disk whose 
density falls exponentially with increasing height (z) 
above the plane as e"^/'*"""''. Its scale height hf^^'^^ ~ 
I ± 0.3 kpc dReid fc Majewski 1993| ; pjha et al. 1996 ; 



Buscr et al. 1999) is about three times larger than that 
of the much brighter thin disk. 

The scale lengt h of the thick disk is h''^'''^ ~ 3±1.5 kpc 
(Buser et al. 1999), similar to that of the Galactic thin 



disk. Despite this, the thick disk contibutes only 2-9% of 



the total local stellar disk light ( 


ieid & Majewski 1993; 


Ojha et al. 1996; 


Buser et al. 199^ 


), and perhaps ~13% 



of the total disk luminosity of the Milky Way 



( [Morrison et al. 19941 ). 

For external galaxies, morphology determined through 
integrated surface brightness photometry is the only cur- 
rent method to detect and characterize faint galactic 
components. Detections of extended light that are 
perhaps indicative of a thick disk component with 
^thick ~ 1 — 2 kpc have been reported in a few external 
edge-on galaxies. Early detections of extra-planar light 
in excess of that as sociated with a thin exponential disk 
were limited to SO ( Burstein 1979) and early- type spirals 
with significant bulges (Ivan der Kruit & Searle 1981a; 



van der Kruit fc Searle 1981b| ; 



Wakamatsu fc Hamabc 1984; Bahcall fc Kylafis 1985 



Morrison et al. 19971) 



Shaw & Gilmore 1989|; |de Grijs fc van der Kruit 19961 ), 

leading to the suppo- 
disks were found in older 
significant central con- 



sition that thick 
stellar systems with 



centrations 



(van der Kruit fc Searle 1981a 



Hamabe fc Wakamatsu 1989, de Grijs fc Peleticr 1997) 



This hypothesis is consistent with the lack of a thick 
luminous component around the small, Scd spiral 
NGC 4244 in deep i?-band observations reaching to 
R — 27.5mag/sq arcsec ( Fry et al. 1999]), and in th e 
bulgeless Sd edge-on NGC 7321 ( [Matthews et al. 19991 ). 
On the other hand, observations indicate that there 
are individual exceptions. Mult iband photometry of the 
later-type Sc spiral NGC 6504 ( |van Dokkum et al. 1991) 
revealed extended light interpreted as a weak thick disk 
with /i^^'^^k ^ 2 kpc. 

Faint light high above the plane of the well-studied, 
late-type, edge-on spiral NGC 5907 has further compli- 
cated the picture of extra-planar light in small- or no-bulge 
spirals. First detected at heights of 3 to 6 kpc above the 
plane in deep i?-band observations (Morrison et al. 1994), 
this extended emission is intriguing because it is unlike 
any known thick disk or stellar component, having instead 
a morphology similar to that inferred for the dark mat- 



ter halo distribution of NGC 5907 ( [Sackett et al. 1994] ). 
Other workers have confirmed the presence of the faint 
light in other bands (BVRIJK), and showed that the 



extended light is redder than the thin stellar disk. If 
the faint light is due to a thick disk, it is unlike any 
other, with a scale length that is at least twice that 



of its thin disk (Morrison 1999). The stellar population 



responsible for this faint light remains highly contro- 
versial, ranging from normal or metal-rich populations 



with steep IMFs ( 


Lcqucux et al. 1996; Rudy et al. 1997 


pames fc Casali 1998 


), old, metal-rich accreted popula- 



tions with normal IMF (Lequeux et al. 1998), or exceed- 



ingly metal-poor or giant-poor populations with few re- 



solvable stars at the tip of the RGB (Zepf et al. 2000). 
The controversy remains because the full spectral energy 
distribution is apparently inconsistent with any s ingle ex- 



planation (e.g. IZepf ct al. 20001 ; |Yost et al. 2000D 



The puzzling nature of the extended light in NGC 5907 
has motivated new studies to test a possible connec- 
tion between faint optical and IR light and dark mat- 



ter in this and other spirals (Gilmore fc Unavanc 1998 



Rauschcr ct al. 1998 


;|Ucmizu et al. 1998 


; Abe et al. 1999 


Bcichman et al. 1999 


Yost et al. 200C; 


Zepf et al. 2000) 



The optical results are mixed, but infrared surface bright- 
ness photometry indicates that whatever produces the 
faint optical light detected to date does not appear to emit 
strongly at IR wavelengths far from the plane of the galac- 
tic disks. Thus, if associated with known stellar popula- 
tions, the sources of the faint light are unlikely to account 
for the dark mass of spiral galaxies. 

In this paper, we report on the collection, reduction 
and analysis of ultra-deep surface photometry of the iso- 
lated, edge-on, low surface brightness, Sd galaxy ESO 342- 
G017, using some of the first science observations taken 
with the VLT. The simple optics, good seeing, and ex- 
tremely well-sampled PSF of our observations ensured a 
low and well-understood level of scattered light and accu- 
rate identification of contaminating sources. Concurrent 
deep observations of unrelated blank fields with the VLT 
were used to create dark sky flat fields at the appropri- 
ate wavelengths. Considering all sources of uncertainty, 
including those from light scattered through the wings of 
the PSF, we conclude that the resulting surface photom- 
etry is reliable to a level of i? = 28.5mag/sq arcsec and 
V — 27.5mag/sq arcsec. Analysis of these data reveals a 



^ The discovery of a faint, long, ver y narrow arc of l ight 
apparently associated with NGC 5907 (Shang et al. 1998) led 



Zheng et al. (1999) to suggest that the extended light in the 
galaxy might be an artifact due to confusion from the arc and 
foreground objects. The arc clearly contributes some light to 
some positions near the galaxy, but is too narrow and asym- 
metric to be the cause of the symmetric extended light de- 
tected by Morrison et al. (1994). Zheng et al. (1999) report 
that their photometry suffers from systematics at light levels 
fainter than R — 27mag/sq arcsec. (Due to a large pixel size, 
the PSF was often undersampled, despite the seeing of 3.4 to 
5.4" that was typical of their observations). Since all detections 
of faint extended light in NGC 5907 have been reported for 
R > 2 7mag/sq arcsec, a nd all optical photometry (including 
that of IZheng et al. 1999) agrees above this level, the mystery 



of this faint halo light remains. 
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Parameter 



Value 



Reference 



ct,5 (J2000.0) 
type 
redshift 
inclination 
PA 

major- axis D* 
ms 
xav 
ran 
m/ 

Mr 

Mv 



21 12 10.8, -37 37 38 
Sc+6 
7680± fOkms-i 
88° 
120? 4 ± 0?5 

86" 
16.67±0.09 
16.40±0.03 
15.92±0.04 
15.47±0.06 
-19.1±0.3 
-18.7±0.3 



Carachentscv ct al. 1999| 
/lathcwson Ford 19961 
/lathewson fc Ford 1996| 
this paper 
this paper 
this paper 



Lauberts fc Valentijn 1989| 



this paper 
this paper 



Mathewson fc Ford 1996 



this paper 
this paper 



*major-axis diameter measured from the R=27.0 mag/sq arcsec contour. 
Magnitudes are not corrected for extinction. 



faint component that we interpret as a thick disk, to our 
knowledge the first thick disk discovered in an LSB galaxy. 

In Sect. 2 we describe the VLT observations and ob- 
serving strategy. In Sect. 3 the data reduction process, 
including the production of dark sky flats and the proce- 
dures for masking, mosaicing, calibrating, and determin- 
ing the sky flux are outlined. The procedure to extract 
profiles from the deep images is given in Sect. 4, along 
with a brief description of the error analysis, which is dis- 
cussed in depth in the appendix. The resulting V and 
R surface photometry of ESO 342-G017 are presented in 
Sect. 5, along with a description of the fitting procedure 
for the thin and thick disk parameters. A thorough anal- 
ysis of scattered light due to the tightly-constrained PSF 
is discussed in Sect. 5, and ruled out as the cause of the 
faint extended hght we detect in ESO 342-G017. The thin 
and thick disks, including their inferred intrinsic proper- 
ties are described in Sect. 6. We summarize and conclude 
in Sect. 7. Throughout this paper we assume a distance 
of 102 Mpc to ESO 342-G017 (based on a Hubble con- 
stant of IIo=75 km/s/Mpc), which yields an image scale 
of 0.495 kpc per arcsecond. 



2. Observations 

2.1. The Target Galaxy 

The target, ESO 342-G017, is a nearby, edge-on galaxy, 
selected on the basis of its right ascension and declination, 
very high disk inclination, absence of a prominent bulge, 
low extinction correction, and optimal angular size. The 
latter is important in order to adequately resolve the disk 
scale height while maintaining sensitivity to faint surface 
brightness in the halo. Our deep i?-band image obtained 
with the VLT-UTl test camera is shown in Fig. nland the 
basic properties of the source are given in Table^. 




-60 -40 -20 20 40 
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Fig. 1. Contour plot of ESO 342-G017 showing levels 
from 20.0 to 27.0 R mag/sq arcsec in 0.5 mag/sq arc- 
sec steps. The image is a central subsection of our total 
ESO 342-G017 mosaic. The inability to trace smooth con- 
tours at the lowest light levels and the noisier background 
on the western end of the source is due to fewer frames 
making up the mosaic on this side of ESO 342-G017. 



2.2. Observing Strategy and Resulting Data 

Bessel V and R observations of ESO 342-G017 were 
made on the nights of 18, 22 and 25 August 1998 as 
part of the ESO VLT-UTl Science Verification (SV) 
program. A complete description of the VLT SV pro- 
gram telescope and instrument set-up can be found in 
Leibundgut et al. 1998 and Giacconi et al. 1999. We give 
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Table 2. Summary of Observations 



Field 


Filter 


Dates 


Total Integration 


Median Seeing 







August 1998 


(seconds) 


(arcsec) 




Bessel- v 


ZZ 


Q Qfin 
ooUU 


-1// 1 
i. i 




Bessel-J? 


1 o oo o cc 
iO, ZZ, Zo 


1 noon 
iUo2U 


u. y 


Flat-field frames: 










HDF-S° 


Bessel-y 


18, 22, 23, 26, 27 


16200 


o'.'g 




Bessel-J? 


18, 22, 23, 25, 26 


15300 


I'.'O 


EIS0046-2930 


Bessel-V 


17 


2700 


O'.'S 




Bessel-J? 


17 


2700 


C.'s 


EIS0046-2951 


Bessel-V 


22 


2700 





only a summary of the issues important for our observa- 
tions of ESO 342-G017. 

The VLT Test Camera, an engineering grade Tektronix 
2048^ CCD, was rebinned 2x2 to improve its surface 
brightness sensitivity, resulting in a binned scale of 0.091 
arcsec pixel"^ and a field of view of 93 arcsec on a side. 
The camera was rotated approximately 60 degrees in or- 
der to position the galaxy major axis along the x-axis of 
the detector. For economy of prose throughout the pa- 
per, we will refer to the northeast and southwest sides of 
ESO 342-G017 as the "northern" and "southern" sides, 
respectively. 

A challenge to our data reduction was the fact that the 
Test Camera CCD is not a science-grade device. As such, 
it displays more than the customary number of cosmetic 
flaws, most noticeably, a large region (^120x130 pixels) 
near the center of the chip with a lower sensitivity than 
its surroundings. Although this "stain" has a strong colour 
dependence (it is more prominent in the blue), we found it 
to be temporally stable and therefore easily corrected with 



our science frame flat- fields (see Sect. 3.2). Furthermore 



ESO 342-G017 was always positioned well away from this 
feature. 

Fortuitously, long total integrations were made of the 
HDF-south and two EIS cluster candidate fields on the 
nights of 17, 18, 22, 23, 25, and 26 August, in the same 
filters as our observations. Using these images to create 
our deep sky flat-fields obviated the time-intensive strat- 
egy of observing off-source fields for ESO 342-G017. Each 
of the images used to make our superflats, as well as the 
observations of ESO 342-G017 itself, were dithered on av- 
erage by more than 10" in both a and S. This allowed for 
the removal of cosmic rays from our galaxy field, and the 
removal of stars in the super sky flat (see Sect. 3.2). 



3. Data Reduction 

3.1. Bias and Dark Current 

The basic image reduction was done using MIDAS. 



The bias frames showed a fixed structure with an over- 
all level that varied up to 20 counts during the course of 
each night. We corrected for this by using the overscan 
region of the detector, which mirrored the same variation. 
For each night, a median- filtered master bias was made 
from at least 20 individual bias images. An average bias 
level was determined for each image from its overscan re- 
gion. The associated master bias was then scaled to each 
overscan mean and subtracted from each image, with the 
0.5 count difference between the overscan and the bias av- 
erage taken into account. No significant dark current was 
measured in the VLT test camera. 



3.2. Creating the Super Sky Flats 

The greatest potential source of error in our final images is 
uncertainty in the flat-field. As many sky counts per pixel 
as possible are required to reduce the statistical error in 
the flat-field which, to avoid large systematic uncertain- 
ties, should be obtained using light with the same spectral 
energy distribution as the primary observation. This was 
done by creating a super flat-field made from careful com- 
binations of the deep EIS and the HDF-S fields that were 
interleaved temporally with our observations of ESO 342- 
G017. The advantage of this method lies in the large total 
exposure of these deep fields, which are devoid of bright 
stars and were well-dithered between individual exposures. 
The HDF-S and EIS fields are located 26?3 and 53?8 away 
from ESO 342-G017, respectively. 

Each candidate sky flat image was inspected visually; 
only those free of defects and temporally close to our ob- 
servations of ESO 342-G017 were chosen. Observations 
of the HDF-S made on 28, 29, and 31 August 1998 were 
not used in our flat-field due to increasing sky levels from 
a waxing moon. The remaining 26 i?-band and 31 V- 
band flat frames contained a total of 73560 and 39550 sky 
electrons per pixel, respectively. Considering only Poisson 
statistics of sky electrons, the flat-field formed from these 
frames should contribute a pixcl-to-pixel error of 0.37% 
(i?-band) and 0.50% (y-band). Of course, variations in 
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the sky brightness across the image and remnant halos 
from inadequately removed bright stars, create large-scale 
errors above that expected from simple Poisson variations. 
We empirically determine the size of this dominant flat- 
field error below. 

The super flat-field was created for each filter sepa- 
rately as follows. Each individual fiat-field sky frame was 
normalised to its modal value as determined in the cen- 
tral 3/5 of the image. The average value of pixel (i, j) was 
then determined from the stack of sky frames for the filter, 
accepting a pixel (i, j, k) from the fcth frame in the compu- 
tation of the average only if it passed two tests. First, its 
deviation from the mean pixel value in the stack at 
must not exceed a given threshold measured in units of 
the noise at that pixel position (a n-a clip) . This criterion 
effectively removed cosmic-ray events and, since each im- 
age was dithered by at least 10" (110 pixels) in both a and 
5 between successive exposures, the bright cores of stars 
and galaxies as well. Second, a median-filtered frame was 
created over a 3x3 pixel window from the average frame 
resulting from the first step. A n-a clip was again applied 
to each pixel {i,j,k) based on the value of its local me- 
dian. The second test was applied to remove any remnant 
faint extended wings of stars and galaxies, which would 
otherwise contaminate the resulting flat-field frame. Only 
pixels satisfying both these "filters" entered the average 
for the flat-field frames. A normalization level was calcu- 
lated from the median value in the central 3/5 of each 
flat-fleld frame, and each image was then flattened and 
renormalizcd. 

In order to test the quality of the flat-fields, and to 
compute an empirical large scale fiat-field error, we re- 
peated the above procedure using only one-half of the 
available HDF-S and EIS images. In this way, flatRl was 
made from HDF-S and EIS images from nights 17, 18, 
22, and 23 August, while flatR2 was made from HDF-S 
and EIS images from nights 23 and 26 August. The two 
subflats Rl and R2 have approximately the same flux lev- 
els. Two y-band subflats were created in the same way. 
The flat-fleld frames flatRl and flatVl were then flat- 
tened using flatR2 and flatV2, respectively. Each was 
then examined visually for any remnant features, and 
then rebinned to a number of relevant scales and the 
rms variation across the frames measured. The cosmetic 
flaws inherent in the Test Camera CCD, particularly the 
"stain" mentioned in Sect. |2.2| , were removed effectively 
by our flat-fleld procedure. The results are summarized 
in Table ^ in which the measured rms is compared to 
that expected from photon statistics alone. The empirical 
values are used in our computation of flat-fielding errors. 

3.3. Mosaicing and Masking the Galaxy Frames 

A region of sky 2f8x2.'2 (R) and 2f2x2f0 {V) around 
ESO 342-G017 was tiled with VLT test camera exposures 
and then combined into a final mosaic. Centroids of a num- 
ber of stars and galaxies (usually 6 to 10) were measured in 
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Fig. 2. Histogram of the number of objects in our i?-band 
image detected by the SeXtractor program, as a function 
of the object's classification. The dividing line between a 
stellar and an extended detection is approximately 0.8; 
the VLT field surrounding ESO 342-G017 is clearly dom- 
inated by background galaxies. 



each individual image to compute their positional offsets 
within the mosaic. In order to remove cosmic ray events, 
images were divided into groups of four closely overlapping 
frames. Using the computed offsets, each group was com- 
bined into a temporary median-filtered image. Each input 
images was compared to its group median and all pixels 
deviating by more than 3.5f7 were replaced by the me- 
dian value. Since cosmic ray events are often surrounded 
by lower brightness halos or tails, a second iteration was 
done at each position at which a cosmic ray was detected. 
In this second pass, a lower pixel correction criteria of 2. Oct 
was applied. 

The 14 (i?-band) and 11 (V^-band) frames with the 
best seeing were then combined, using integer pixel shifts, 
into R- and V-band mosaic frames. Given the small pixel 
size and large over-sampling, this did not limit the resolu- 
tion of our resulting image. Since different regions of the 
mosaic are constructed from different numbers of images, 
it is necessary to renormalize. To do this an identical set 
of frames was created having the same sizes and offsets, 
but containing only the modal value of the source-free sky 
background. These were also combined into a mosaic and 
used to renormalize the R- and l^-band mosaic frames. A 
subsection of the i?-band image resulting from this proce- 
dure is shown in Fig. |l|. 

In order to be able detect faint light associated with 
ESO 342-G017 in our deep mosaic, foreground stars and 
background galaxies must be masked out. Since ESO 342- 
G017 was explicitly chosen for its paucity of foreground 
stars, most of the objects contaminating its background 
are galaxies (see Fig. ||) , and simple profile fitting cannot 
be used to model and subtract contaminants. 

Instead, we used the SeXtractor detection algorithm 



(Bertin & Arnouts 1996) to find sources not associated 



6 Neeser et al.: Thick Disk in LSB Spiral ESO 342-G017 



Table 3. RMS Flatness of Flat-fields 



Flat Correcting Flat 


Filter 


Rebinned Size (") 


Relevant Scale 


Measured rms 


Pixel — to — Pixel rms 


flatRl/flatR2 
flatVl/flatV2 


R 
V 


0'.'091 
0'.'091 


1 pixel 
1 pixel 


0.57% 
0.78% 




tlatRl/tlatR2 
flatVl/flatV2 


R 

V 


0'.'806 

o'.'soe 


400 pc (~hdisk) 
400 pc (-hdisk) 


0.11% 
0.14% 


0.064% 
0.088% 


flatRl/flatR2 
flatVl/flatV2 


R 

V 


o'.'g 
I'.'i 


450 pc (PSF FWHM in R) 
550 pc (PSF FWHM in V) 


0.16% 
0.12% 


0.058% 
0.065% 


flatRl/flatR2 
flatVl/flatV2 


R 
V 


6'.' 04 
6'.' 04 


3kpC (-^hhalo) 
3kpC (~hhalo) 


0.08% 
0.11% 


0.0086% 
0.012% 




Fig. 3. Final i?-band (top) and y-band masked images 
of ESO 342~G017. Objects detected with SeXtractor in 
either band have been masked in both frames. Levels 
— 3.5asfey to lOCTsfey around the frame median are shown, 
where Osky is the background rms/pixel of the frame. 

with ESO 342-G017. A source was defined to consist of at 
least five connected pixels at a level of l.ha above the lo- 
cal background, which was computed over a 32x32 pixel 
mesh. The so-called OBJECTS output of SeXtractor, es- 
sentially a frame of all detected objects separated by null 
pixels, proved valuable in creating a mask for objects be- 



yond the outermost contours of ESO 342-G017. The initial 
output masks still retained a faint halo of emission around 
brighter sources. For this reason, the masks were grown in 
size iteratively until a histogram of the unmasked back- 
ground pixels no longer changed shape, indicating that the 
local background level had been reached. 

A crucial step in the data reduction process is the de- 
termination of an accurate value for the background sky 
level. A large central section of both masked mosaics was 
extracted so that its area contained the largest possible 
number of overlapping individual images (>11 for the R- 
band and >8 for the V-band). In order to prevent any 
emission from ESO 342-G017 contributing to the sky sig- 
nal, the galaxy was liberally masked out to 20'.' 1 (lOkpc) 
above and below the central plane of its disk, and along 
its major axis to the outermost edges of the images. The 
mask sizes of the brightest field stars were also liberally 
increased for this procedure. 

3.4. Determination of Sky Level 

The distribution of sky values are shown in the histograms 
of Fig. 0, which were used to compute the true background 
value of the unmasked pixels in each image, and the as- 
sociated error in its mean. These sky values are Sr — 
16651.5 ±0.4 e- pix^^ and Sy = 2950.2 ±0.2 e" pix"! 
in the R and V bands, respectively. Using the calibra- 
tion described in the next section, these values corre- 
spond to msky(R) — 20.98 mag/sq arcsec and msky(V) — 
21.60 mag/sq arcsec, with a systematic uncertainty domi- 
nated by calibration errors of ~5%. The systematic devia- 
tion from gaussian behaviour seen at extreme pixel values 
in Fig. 1^ is slight and very much smaller, in its integrated 
effect on the average sky value, than the uncertainties 5S 
based on gaussian statistics reported above. 

3.5. Calibration to Standard System 

Our photometric calibration was based on results supplied 
by the SV team together with the distribution of our data. 
A photometric solution was available only for the observa- 
tions of ESO 342-G017 on 22 and 15 August, as these were 
the two photometric nights. Typically, four standard fields 
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Fig. 5. The positions of the profile extractions shown on the mosaiced, masked, i?-band image. The F-band image 
was extracted at the identical positions, but since it is smaller (see Fig. ||), the V profiles only reach number 52. The 
vertical profiles averaged together to create Figs. and ^ are labelled at the top. 



were observed several times during each of these nights, 
with an average of about 10 Landolt standard stars being 
used to compute the photometric solutions. The standards 
chosen spanned a significant range of colours in order to 
adequately measure the colour term. 

4. Extracting Vertical Profiles 

Achieving acceptable signal-to-noise at surface brightness 
levels 6 to 8 mag/sq arcsec below sky requires averaging 
over a large number of pixels. We begin by extracting a 
number of vertical rectangular regions, each of dimension 
21 X 530 pixels (0.9x 24kpc), perpendicular to the disk 
of ESO 342-G017. These extracted areas are centered on 
the major axis of the galaxy, avoid the most prominent 
HII regions, and extend well beyond the visible disk. The 
positions of the extractions were identical for both the R 
and y-band images; 71 regions were extracted from the 
i?-band image and 52 from the smaller y-band image. 
Fig. H shows these areas atop on our masked mosaic R- 
band image. 

From these initial extractions four levels of averaging 
were performed in order to increase the signal-to-noise: 

1. The sum of the flux across the 21 pixel wide x-direction 
was determined for each extraction, and normalized 
by the number of non-masked pixels contained in each 
row. 

2. Due to the extremely symmetrical cross-section of 
ESO 342-G017 (see Fig. ^) we were able to average 
the profiles above and below the disk. 



3. In the vertical direction (z-axis), each profile was av- 
eraged over the size of the seeing FWHM (~10 pixels) 
for all points above the plane of the galaxy. 

4. Finally, a number of profiles were averaged together. 
An average of the three innermost profiles (extractions 
36 to 38) was made to create one central profile, groups 
of four profiles (28-31 and 32-35 to the east and 39- 
42 and 43-46 to the west) were averaged on each side 
of the center, groups of six (16-21 and 22-27 to the 
east and 47-52 and 53-58 to the west) were averaged 
on the outermost ends of the galaxy. The extractions 
averaged together are shown at the top of Fig. ||. In 
the case of the V-band image, the averaging process 
ends with extraction number 52. 

The resulting masks made from the R and V images 
separately, were then multiplied together to create a mas- 
ter mask frame that was applied to each mosaic. This pro- 
cedure masked 10.3% and 11.0% of the total image areas 
in the final R and V mosaics respectively. The masked 
images are shown in Fig. |3| 

For each of the vertical extractions covering the visible 
disk of ESO 342-G017 (profiles 16 to 58 in R and 16 to 
52 in V), a least-squares fit to the thin disk component 
was made. A simultaneous two-component (thin and thick 
disk) fit was made to each extraction, using an exponential 
parametrization given by 

f{z) = C'^exp{-\z\lhfn + C''^^exp{-\z\lhf'^'^) (1) 

for both components, where /□ is the surface flux at 
the position the extraction crosses the major axis of the 
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Table 4. Typical Errors for ESO 342-G017 Vertical Profiles 



Uncertainty 


R-band 




V-band 




[electrons] 


Galaxy Center (%) 


±4kpc (%) 


Galaxy Center (%) 


±4kpc (%) 


Averaged flux per pixel 


30640 


16703 


5498 


2958 


Sky flux per pixel 


iODOi.O 


10001.0 


zyou.z 


zyou.z 


Net flux per pixel 


13988.5 


51.5 


2547.8 


7.8 


Read Noise (ctrn) 


0.04 (0.0003) 


0.04 (0.08) 


0.05 (0.002) 


0.05 (0.6) 


Flat-Fielding (ctff) 


7.4 (0.05) 


4.0 (7.8) 


1.8 (0.07) 


1.0 (12.5) 


Photon Noise (ctpn) 


1.0 (0.007) 


0.8 (1.5) 


0.5 (0.02) 


0.4 (5.4) 


Mosaicing Error ((Jm) 


1.7 (0.01) 


0.9 (1.8) 


2.6 (0.1) 


1.4 (17.9) 


Surface Brightness Fluctuations (ctl) 


0.2 (0.001) 


0.01 (0.02) 


0.06 (0.002) 


0.003 (0.04) 


Total Statistical Error (ctstat) 


7.6 (0.05) 


4.2 (8.1) 


3.2 (0.1) 


1.7 (22.3) 


Sky Subtraction (ass)* 


0.4 (0.003) 


0.4 (0.8) 


0.2 (0.008) 


0.2 (2.6) 


Total m ± Am [mag/sq arcsec] 


21.17 ± 0.04 


97 o +0.13 
^ ' -0.12 


21.76 ± 0.03 


28.0 



The numbers in parentheses are the errors as a percentage of the sky-subtracted flux. 
* Note: the sky subtraction error is a systematic error and, therefore, does not enter into ctstat- 




Fig. 4. Histogram of sky pixels from the completely 
masked i?-band (top panel) and F-band (bottom panel). 



A gaussian fit to the histograms {exp{—{x — xq) /2a^) 

a{R)=120.2 e~i pix^^ 



with xo(-R)=16651.5 e" 
and xo(F)=2950.2 e"! 
shown as a dotted line. 



pix 
pix"\ 



cr(t/)=39.7 e-i pix- 



is 




galaxy, z is the projected distance from the major axis, 
and hz is the exponential scale height. The fitted param- 



IDistance above/below plane [kpc] 

Fig. 6. The symmetry of the vertical surface brightness 
profiles of ESO 342-G017. The north (dashed fines) and 
south (dotted lines) profiles extracted at various positions 
along the disk. The object-masked images have been used. 

eters are hf''\ f^^'"^, and hf'"^. Regions strongly 

affected by dust or clumpy HII regions were excluded from 
the fit. For comparison, a single-component (thin disk) fit 
was also performed for those data that lie between 1 and 
3 mag/sq arcsec below the central galaxy surface bright- 
ness. Results are presented in Sect. ^. 

4.1. Error Analysis 

We present here a brief summary of the sources of pho- 
tometric uncertainty and their magnitudes; the reader is 
referred to the Appendix for a more detailed discussion. 

For illustration, Table ^ shows the average fiux lev- 
els and uncertainties at two positions along the central 
vertical profile (the average of extractions 36 to 38) of 
ESO 342-G017, one at the galaxy center and another at 
the much fainter light levels 4 kpc above the galaxy disk. 
For each of the two flux extremes in each of R and V- 
bands, we give the uncertainties associated with the av- 
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erage flux per pixel (averaged over the unmasked area of 
3x21x 10 pixels) in units of electrons and as a percentage 
of the sky-subtracted flux (given in parentheses) . 

The systematic uncertainties in the sky level of SSr = 
0.4 pix~^ and SSy — 0.2 pix^^ correspond to 
errors of only 0.0024% and 0.0068% per pixel in the R 
and V bands, respectively. These systematic uncertainties 
are present in the sky-subtracted proflles we present in 
the next section, but because they correspond to light 
levels AR = 11.5 and AV = 10.4 magnitudes below 
the sky (msky(R) = 20.98 mag/sq arcsec and msky(V) ~ 
21.60 mag/sq arcsec), they are of no importance (<10%) 
over the range of surface brightness we consider. The sys- 
tematic uncertainty in overall calibration to a standard 
system of about 5% is relevant, but simply corresponds to 
a possible overall shift in the surface brightness scale by 
that amount. Note that at bright flux levels, the error in 
the magnitudes is dominated by the error in the photo- 
metric conversion term, not by crgxAT- At faint flux levels, 
the situation is reversed. 

5. Results 

5.1. Analysis of the Surface Brightness Profiles 

In Figs. ^ and ||, we show the vertical extractions, av- 
eraged above and below the plane of ESO 342-G017, de- 
rived from our deep VLT imaging and discussed in Sect. ^. 
Individual extractions consist of the horizontal average of 
21 pixel wide rectangles, with foreground stars and back- 
ground galaxies masked. Extractions above and below the 
galaxy disk were averaged about their axis of symmetry 
and averaged in groups to produce the profiles shown in 
the figures. In order to display meaningfully our data at 
the faintest levels, insets in Figs. and || show fluxes on 
a linear scale for distances greater than 6 kpc from the 
major axis of ESO 342-G017. The scatter in each inset 
about zero indicates clearly that the sky flux has been 
well-subtracted in our flnal mosaic within our calculated 
uncertainties. 

The deviation from pure exponential behaviour in 
nearly all of the profiles indicates that presence of ex- 
tended light beyond that expected for a purely exponen- 
tial stellar disk. This motivated our choice of one- (thin 
disk only) and two-component (thin-|-thick disks) least- 
squares fits to the profiles, which are overplotted in Figs. 
and H on the data. Both components were modeled as ex- 
ponential disks (see equation 1) with the scale height 
and central surface brightness /o as free parameters. The 
scale heights and central surface brightnesses (expressed 
in mag/sq arcsec) derived from the simultaneous thin plus 
thick disk fits are summarized in Table ^ 

5.2. Could the Faint Extended Emission be Scattered 
Light? 

The first-order effect of turbulence in the atmosphere 
causes the radial point spread functions (PSFs) of point- 



like objects measured by an astronomical detector to have 
a roughly gaussian shape, but many effects, including 
scattering in the telescope optics, can lead to broader 
wings. Although the simple optics of the VLT test camera 
(Giacconi ct al. 1999) should minimize such a scattering, 
faint wings in the PSF are still present. To quantify the 
effect of these wings on our faint surface brightness pho- 
tometry, we measured the PSFs of isolated fainter stars 
in the held of ESO 342-G017 and bright standard stars 
observed on the same nights as our science frames. 

In order to be meaningful, such PSFs must be con- 
structed with high signal-to-noise data. A conservative es- 
timate of the precision required can be made by assuming 
that aU of the hght of ESO 342-G017 is conflned to a point 
at a distance equal to the angular separation between the 
center of the galaxy and the most distant point above the 
plane we consider. Such an estimate shows that ,at 6 kpc 
above the plane, the amount of i?-band light brighter than 
30 mag/sq arcsec scattered from ESO 342-G017 can be 
quantifled easily if the PSF is known to a precision of 
~ 2.5 X 10~^ at that distance. 

Only three relatively isolated stars near the center of 
the ESO 342-G017 mosaic are available; the R- and V- 
band images of these were masked, added, and azimuthally 
averaged. The result is shown in Fig. ||. Note that although 
the extended emission around ESO 342-G017 is seen more 
clearly in the i?-band, median seeing in R is better than in 
V. At faint light levels, both the R and V faint-star PSFs 
have broader wings than a pure gaussian. Unfortunately, 
the statistical noise in these faint-star PSFs, and the rela- 
tive size of the systematic photometric uncertainties over 
the relevant radii, precludes measurement in the wings to 
the accuracy we require. In principle, saturated stars on 
the mosaic could be used to study the wings of the PSF, 
but our small fleld contained only two; one has a near 
bright neighbor and the other does not fall on the y-band 
mosaic. 

We therefore study the PSF wings using much brighter 
standard stars imaged during the same observing run. We 
build a model PSF directly from the data, using the iso- 
lated three stars on the mosaic of ESO 342~G017 to derive 
the PSF out to 2.5", and a bright reference star observed 
with similar seeing to derive the PSF from 2.5 to 16" (i.e., 
out to 8 kpc above the galaxy plane) . Since the seeing was 
slightly better during the imaging of the reference star, the 
standard star profllc was horizontally displaced to create 
a smooth match to the inner PSF derived from the mo- 
saic. The result is shown as the thin solid line in Fig. |[ 
The i?-band PSF of the standard star is consistent with 
zero at the level of 2.5 x 10^^ from 7 to 12 arcsec {^-^ 3.5 
to 6 kpc), satisfying the conservative requirement that we 
derived above. We conclude, therefore, that the PSF is 
well enough understood to determine its effect on the ob- 
served shape of the vertical surface brightness proflles of 
ESO 342-G017. 

In order to examine whether the extended light ap- 
parent in Figs. 1^ and || might be due to thin disk light 
scattered through the broad wings of the PSF to other 
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Fig. 7. The i?-band averaged profiles through the disk of ESO 342-G017, perpendicular to the major axis. The average 
position from the galaxy center is given at the top of each panel (east of center is indicated by i? >0; west of center by 
R <0). The insets show, on a linear scale, the background-subtracted flux levels of each profile at distances more than 
6 kpc from the galaxy disk. The simultaneous thin and thick disk fit and the range of data used for the fit is shown 
by a solid line; the dotted line is the extrapolation of the fit. The dashed line is intended as a guide, and indicates 
a single-component fit to data dominated by the thin disk. This fit was restricted to data between 1 and 3 mag/sqr 
arcsec fainter than the peak flux. 



positions on the detector, we convolved a model exponen- 
tial disk with intrinsic structural characteristics similar to 
those of ESO 342-G017 with our model i?-band PSF. The 
intrinsic thin disk model parameters reported in Sect. ^ 
were determined by requiring that, after inclination and 
convolution with the observed PSF, the projected thin- 



disk fitted parameters were retrieved. The degree to which 
the thin disk fits are reproduced is illustrated in Fig. |o|. 

Due to its high inclination, ESO 342-G017 has an ob- 
served surface brightness along its length that is much 
larger than the intrinsic (input) face-on value. Except for 
the central regions, which suffer a net loss of light from 
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Fig. 8. As in Fig. 0, but for the ^-band averaged profiles through the disk of ESO 342-G017. 



scattering, the primary effect of incHnation — and to a 
lesser extent scattering — is to increase the amount of 
light observed at a given angular distance from the plane 
of ESO 342-G017. The result (output) is an observed pro- 
file that is approximately exponential, but with a pro- 
jected scale height larger than the intrinsic value. More 
importantly, however. Fig. |l^ clearly illustrates that for 
surface brightnesses brighter than R « 28.5 mag/sq arc- 
sec, no substantial light in excess of the projected thin 
disk profile is generated by inclination and scattered light 
effects. The extended light R > 26.5 mag/sq arcsec in 
many of the profiles of Fig. therefore, must have an- 



other cause; we conclude that it is intrinsic to the galaxy 
itself. This conclusion is supported by the constant color 
(or possible slight reddening) of the extended light with in- 
creasing distance from the galaxy plane, despite the fact 
that the scattering in the V^-band images is larger than 
that in R as measured from the stellar PSF on the science 



6. The Thin and Thick Disks of ESO 342-G017 

The extended hght in ESO 342-G017 is reasonably well fit 
by a thick exponential disk with nearly constant projected 
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Table 5. Fitted Disk Structure Parameters for ESO 342-G017 



fj.{0) ( maj 
R (kpc) R 



-17.5 
-11.8 
-7.1 
-3.3 

3.3 
7.1 
11.8 
17.5 



Thin 

;/sq arcsec) 
V 



Disk 



ftz (pc) 



R 



V 



li{0) ( maj 
R 



Thick 
;/sq arcsec) 
V 



Disk 



ftz (pc) 



R 



24.2 
22.3 



1.2 
1.3 
1.2 



21.25 
20.01 



20.83tH3 

+0.12 
-0.10 
-0.04 
-0.04 

20 iq+f 02 
20 61+" °'' 
20.901^6 

+0.51 
-0.34 



7+1.5 
-1.3 

22.lt\i 

+0.09 
0.09 
+0.17 
0.17 
+0.12 
-0.10 
+0.14 
0.12 
0.28 
0.22 



22. 7"^ 



21.37 
20.49 
20 
21.03 
21.05 



21.44 



340±240 
520±190 
421±18 
494±35 
346±9 
396±13 
416±22 
431±33 
324±51 



350±210 
320±210 
410±140 
428±91 
360±46 
384±75 
422±87 
390±110 



22.493 
23.0:! 



(-0.9U 
-0.85 
1+1.9 

23.05l°S 
22.4312:^3 

99 1 1+0.09 
^^•11-0.08 
99 OO+0.22 

99 c;q+0.33 
ZZ.OO_g 25 
-1.1 



22.9 
22.08 
23.0 



1.5 
+ 1.0 
-0.96 
+1.1 
1.1 



23.62 
24.85 



+0.90 



on 7Q + 1.0 

99 c;Q-t-0.57 
z,i.ao_Q 43 

99 9q+i-7 

90 91 +0.90 

-1+1.5 
-1.4 



23.0: 



600± 120 
720± 180 
870± 55 
807± 59 
808± 16 
806± 33 
856± 52 
950± 160 
1100± 220 



V 



600±230 
600±160 
760±270 
820±140 
806±84 
740±120 
880±260 
760±320 



Note: these structure parameters are the results of the model fits to the observed profiles and have not been corrected for 
projection effects. 



scale height (hz) as a function of galactocentric radius (R) 
along the major axis of the galaxy. This is illustrated in 
Figs. and ||, and in Fig. where the fitted values of 
hz for each averaged extraction are shown in both the 
R and V bands for the thick and thin exponential disk 
components. The error-weighted mean of the projected 
scale heights are: /i*/'" = 380 ± 35 pc and hf''^^ = 810 ± 
40 pc in the i?-band and /if'" = 380 ± 45 pc and hf''=^ = 
760±75 pc in the TZ-band. The projected scale length, /i^, 
of the thin disk is more difficult to assess, but is estimated 
from the fitted values of /i(0) as a function of position 
along the major axis to be about 8.9 ± 1.5 kpc in both 
bands. The projected scale length of the fitted thin disk is 
indeterminate from the V-band frames, but is consistent, 
within the uncertainties, with the projected scale length 
of the thin disk in the R-band. 

When deprojected and deconvolved (that is, taking 
into account line-of-sight effects due to the inclination of 
the galaxy and seeing) , the true face-on surface brightness 
of the thin disk in the i?-band is /i^'o ~ ^3.6 mag/sq arc- 
sec, with a true scale height and scale length of h^^p" — 
310 pc and h.^^^ = 5.9 kpc, respectively. These estimates 
were made by convolving model thin exponential disks in- 
clined at 88° with the measured PSF, and requiring that 
the resulting vertical and radial profiles matched those 
fitted to the observed profiles (see fig |l^). The intrinsic 
thin disk scale heights and lengths in the F-band are 
the same as those in the i?-band, within uncertainties. 
The projected scale height, hz, is larger than the intrinsic 
value hz,o primarily due to convolution with the compa- 
rably sized PSF. On the other hand, the projected scale 
length, hfi, is larger than /i_r.o because of line-of-sight ef- 
fects due to the extreme inclination of the galaxy. The 
thin disk has an inferred face-on surface brightness in V 
of 

Mv;o" = 24.1 mag/sq arcsec, implying an intrinsic color 
of V — R « 0.5 for the thin disk. Since the color is found by 
extrapolating the fitted parameters into the plane of the 
disk, it is relatively, though not completely, insensitive to 
dust and clumpy luminosity from HII regions. We estimate 
that the uncertainty in our inferred intrinsic parameters 



is about 10— 15%, primarily coming from uncertainties in 
inclination and the fit parameters. 

The structural parameters of the extended light are 
more uncertain, but also much less affected by inclination 
and seeing effects. We have not attempted, therefore, to 
deproject the thick disk scale parameters, but expect that 
in the R band the intrinsic scale height is close to the pro- 
jected value of hf^^^^ = 810 ± 40 pc, while the true scale 
length of the thick disk is between ^--^6 and 9 kpc, (the 
intrinsic thin disk and projected thick disk values, respec- 
tively). The value of the central surface brightness of the 
thick disk is uncertain, but can be constrained. For a pure 
exponential disk, the edge-on central surface brightness 
Sedge (0) (in linear units) can be shown to be given by 
Sedge(O) = (hR_o/hz,o)So(0), where So(0) is the face-on 
central surface brightness. If we assume that the fitted 
value /i^"^"^ = 22.1 mag/sq arcsec of the thick disk in 
the i?-band is a good approximation to the actual edge- 
on value for the thick disk, then, based on our estimates 
of these quantities and their uncertainties, we can deduce 
that 24.1 < fi*'^^'^^^ < 24.9 mag/sq arcsec. The PSF may 
have a small efifect that would cause the fitted value to 
be higher than the actual value, in which case these con- 
straints would be pushed to slightly fainter magnitudes. 
The detection of the thick disk in the V^-band is less se- 
cure, both because the S/N of our relative surface bright- 
ness photometry is lower in V and because the PSF (and 
thus scattered light problems) is larger in V. Furthermore, 
beyond galactocentric radii of 5 kpc, there is only a small 
statistical difference in the inferred scale heights of the fit- 
ted thin and thick disks (Fig. |ll|), and the extrapolated 
in-plane surface brightness of the thick component in the 
V^-band shows no clear trend with major axis radius. 

The intrinsic i?-band scale heights of the thin and thick 
disk components of ESO 342-G017 are similar to those of 
the Milky Way, but because the intrinsic scale length of its 
thin disk is larger than the commonly accepted Galactic 
value of hjifi = 3—3.5 kpc (see references in Backett 1997), 
the ratio hR,^o/hz,o is --50% larger for ESO 342-G017 than 
for the Galaxy. Since the total luminosity of any pure ex- 
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Fig. 9. Measured radial point spread functions (PSF) for 
standard stars and fainter stars in the ESO 342-G017 field 
shown on magnitude (top) and linear (bottom) scales cho- 
sen to display the full dynamic range in a meaningful way. 
The data shown in the two plots are otherwise identical, 
and the apparent change in PSF shape in the linear plot is 
solely due to the different scales of the logarithmic and lin- 
ear representation. The R and V PSFs taken from isolated 
stars in the ESO 342-G017 field are shown as solid and 
open dots, respectively. The higher signal-to-noise PSF de- 
rived from a brighter standard star is shown as the dashed 
line. The model PSF formed by matching the faint star 
PSF to the wings of the bright star PSF is indicated by 
the solid line. For comparison, a gaussian with FWHM = 
1" (dotted line) is also shown. Error bars shown in the 
lower panel represent photon noise only. 



ponential disk is given by L = 2iTT,Q{0)hj^ g, if the intrinsic 
scale lengths of the thick and thin components are equal, 
the ratio of total light in each is given by the ratio of their 
intrinsic central surface brightness. Together with the con- 
straints on fiRfl for the two components derived above, 
this assumption implies that the thick disk contributes 
~20-40% of the total i?-band hght of ESO 342-G017, ex- 
cluding the light in individual masked HII regions. 

Finally, we note that these constraints on the lumi- 
nosity contribution of the thick disk imply a combined 
(thin+thick disk) face-on central surface brightness for 
ESO 342-G017of^/j,o > 23.1. Since the B-i? color of the 
galaxy is certainly greater than zero, and probably > 0.5, 
this places ESO 342-G017 firmly in the class of low sur- 
face brightness (LSB) galaxies, which are generally defined 
as those disks with B-band face-on central surface bright- 
nesses fiB,o > 23mag/sq arcsec (c.f., de Blok et al. 1995| ). 




and Fitted Profiles 



2 4 6 8 

Distance from Galaxy Center (kpc) 

Fig. 10. The effect of inclination and PSF convolution on 
the observed radial and vertical surface brightness pro- 
files along the major and minor axes of a thin exponential 
disk similar to that of ESO 342-G017. The inferred intrin- 
sic vertical and face-on radial profiles arc shown as thin 
dashed lines. The thick solid line indicates the result after 
inclination by 88 degrees and convolution with the high 
signal-to-noise PSF determined from isolated faint stars 
on the science mosaic and a bright standard star. The 
vertical (minor axis) and radial (major axis) profiles of an 
exponential thin disk with typical fitted parameters for 
the projected scale height, — 380 pc, and scale length, 
hfj = 8.9 kpc, (see Sect. @) are shown as thin solid lines. 



7. Summary and Conclusions 

We have used the VLT test camera on UTl to obtain deep 
surface brightness photometry of the edge-on LSB galaxy 
ESO 342-G017 in the V and i?-bands. Careful masking 
of foreground and background objects to obtain an accu- 
rate value of the sky flux on our science mosaics, and an 
analysis of flat-fielding uncertainties - both statistical and 
systematic - on a variety of spatial scales, allow us to esti- 
mate confidently the total uncertainty in our deep surface 
photometry. We conclude that on the size scales important 
for probing faint, extended structure, we reach ^ = 28 and 
R — 29 mag/sq arcsec. A detailed analysis of the PSF of 
the images, derived from faint isolated stars on the mosaic 
and standard stars, indicates that scattered light affects 
the extended vertical luminosity profiles of ESO 342-G017 
only for R > 28.5 mag/sq arcsec. 

Extended light in excess of that expected for a 
single-component thin disk is detected at about R > 
26.5 mag/sq arcsec in nearly all vertical profiles per- 
pendicular to and up to 17 kpc along the major axis of 
ESO 342-G017. The same component may have also been 
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1200 - • Thin Disk 
o Thick Disk 



10 -10 

R (kpc from disk center) 



1200 - • Thin Disk 
o Thick Disk 



R (kpc from disk center) 

Fig. 11. The fitted values of hz for two-component (thin 

+ thick exponential disks) model of the vertical surface 
brightness extractions of ESO 342-G017 in both R (top) 
and V (bottom) bands. The error bars indicate the formal 
errors of the fit, and arc clearly larger in V-band and at 
larger galactocentric radius R where the S/N is poorest. 
Horizontal dashed lines indicate the error-weighted mean 
of hz for the two components in each photometric band. 



detected in the V band frames, but the lower S/N of these 
frames and the larger PSF in V make this detection less 
robust. Given the geometric form of the extended light 
in this apparently bulgeless galaxy, we interpret the faint 
-R-band light as a thick disk. 

Two-component exponential disk fits were made to the 
observed surface brightness profiles and used to determine 
projected and - after deprojection and deconvolution - in- 
trinsic structure parameters for the thin disk of ESO 342- 
G017 in the V- and i?-bands and for the thick component 
in R. In particular, we find: 



• R band (thin disk) 
oV-band (thin disk) 

* R-band (thick disk) 




R (kpc from disk center) 



Fig. 12. The fitted values of /u(0) for our two-component 
fits to the vertical surface brightness extractions of 
ESO 342-G017. in both R and V bands. The F-band 
thick disk fit is not shown due to its poorer S/N. As a 
rough guide, an exponential of the form f{R) = /(O) * 
exp{^R/hji) has been drawn on the data using a con- 
stant scale length of Hr = 8.9 kpc and the peak central 
surface brightness (^^(0)) of each component. The intent 
is to show that a hji > 8 kpc is consistent with the data, 
and that the Hr is similar for both thin and thick disks. 



— The thin disk has projected scale heights perpendicular 
to the major axis of /i**^'" = 380 ± 35 pc in the _R-band 
and /if'" = 380 ± 45 pc in V. The projected scale 
length of the thin disk is 8.9 ± 1.5 kpc in both bands. 

— After deprojection and deconvolution with the PSF de- 
rived from our observations, we estimate that, within 
the errors, the intrinsic thin disk scale heights arc 80% 
of the measured values. The intrinsic scale length of 
the thin disk is 2/3 of the fitted value. 

— The face-on central surface brightness of the thin disk 
is estimated to be //.q"" = 23.6 mag/sq arcsec in R 
and /^o*^'" = 24.1 mag/sq arcsec in V . 

— The thick disk, which is detected robustly in our R 
surface photometry, has a projected scale height of 
^thick ^ g;^Q ^ 4Q ^-j^g i?-band and hf'''^ = 

760 ± 75 pc in V. The projected scale length of the 
thick disk cannot be determined precisely from our ob- 
servations, but within uncertainties is consistent with 
that of the thin disk. 

— The intrinsic scale parameters of the thick disk are 
somewhat smaller than these measured projected ones, 
though not as dramatically different as the projected 
and intrinsic parameters of the thin disk. 

— Simple considerations lead to an estimate for the face- 
on central surface brightness of the thick disk of 24.1 < 

, ,thicl< 



.thick ^ 24.9 mag/sq arcsec. 
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- The thick disk is hkely to contribute 20-40% of the 
total (old) disk i?-band luminosity of the galaxy. 

— The total central surface brightness of the (thin + 
thick) disk is fiR^o > 23.1 mag/sq arcsec, which places 
ESO 342-G017 securely in the category of low surface 
brightness (LSB) galaxies. 

This detection of a thick disk adds to only a few others 
known in external galaxies (see Sect. 1), and to our knowl- 
edge is the first known thick disk in an LSB galaxy. The 
thick and thin disks of ESO 342-G017 have similar scale 
heights as their corresponding components in the Milky 
Way, but larger scale lengths. Importantly, the thick disk 
of ESO 342-G017 appears to contribute a larger fraction 
of the overall old disk light than does the Galactic thick 
disk. (Young HII regions have been masked and so do not 
enter into the extrapolated estimates we have made.) A 
prominent thick disk in ESO 342-G017 is particularly in- 
teresting since, compared to their high surface brightness 
cousins, LSB galaxies are thought to be more dark-matter 
dominated and to have less evolved disks. The VLT ob- 
servations reported here suggest that, at least in the case 
of ESO 342-G017, such an unevolved thin disk can coex- 
ist with a substantial thick luminous component, perhaps 
providing a clue to the formation mechanism of thick disks 
in all spirals. 
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Appendix A: Photometric Uncertainties 

Accurate detection and characterization of faint surface 
brightness features in galaxies requires a thorough under- 
standing of the uncertainties in CCD photometry at very 
low light levels. Both systematic and statistical uncertain- 
ties are present, and can affect the photometry over differ- 
ent spatial scales. We consider here seven different sources 
of photometric uncertainty, and combine them to create 
an error budget for an area of arbitrary size in the com- 
bined mosaic of the ESO 342-G017 field. The error bars 
presented in Figs.|^ and || are calculated according to this 
error budget. 

Our analysis of photometric uncertainties must reflect 
the process by which the deep, masked mosaics from which 
we derive surface brightness proflles were generated. In 
what follows, all fluxes and uncertainties will be expressed 
as numbers of electrons e~, and i will be used as an index 
to label one of the Nf individual frames that were co- 
added to form the mosaic at that position. 

The flux FsMix,y) at any sky position in the mosaic 
is defined by: 



C 



Nf 



N, 



where Si is the sky value (in e^) on frame i, that is, the 
median of all background (non-object, non-masked) pixels 
in frame i. The quantity is a normalization constant 
defined by 



1 Nj 



I 



N. 



sky 



Nsky 



where j runs over all sky pixels on frame i, and i runs 
over all frames used anywhere in the matrix. On the other 
hand, note that if frame i did not contribute to the flux at 
position {x,y), due to it being masked, then e~(x,y) — 
and Si = 0. Thus Csm is a property of submosaics - por- 
tions of the mosaic composed from the same individual 
frames - and is of the order -/Vg^, where iVsM is the num- 
ber of frames contributing to submosaic SM. The electron 
flux at position (a;, y) is thus the total flux of all frames 
contributing to the mosaic at that position, weighted by 
the total sky value on each part of the submosaic in order 
to account for frame-to-frame differences in transparency 
and total exposure time. 

We now consider separately individual sources of pho- 
tometric uncertainty within detector regions composed of 
Np pixels spread over a total area A, of which are 
unmasked pixels that combine to form an unmasked area 
A'. Throughout, we will assume that all portions of the 
subarea A on the mosaic were constructed from the same 
individual CCD frames. All expressions for uncertainties 
are expressed in terms of numbers of electrons. 
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A.l. Read Noise: grm 



A.3. Photon Noise: apN 



The noise associated with reading the charge collected in 
the CCD array is associated with every pixel of the array. 
For the UTl test camera, this noise has a random distribu- 
tion with an rms (root-mean-square) value of 7.2e~ . Over 
an unmasked area A' on the detector, the uncertainty in 
the flux contributed by read noise is thus 



SM 



NfN; 



assuming that each of the Nf frames contributing to the 
mosiac area A' had Np unmasked pixels. 

A. 2. Flat-Fielding: app 

Flat-fielding was performed by constructing supersky flats 
from moonless UTl Hubble Deep Field South and EIS im- 
ages in the same bands taken during a 10-day period co- 
inciding with our ESO 342-G017 observations. Dithering 
helped to ensure that sky objects did not fall on the same 
portion of the physical detector and could thus be removed 
in the median process (see Sect. 3.1). Nevertheless, Table ^ 
demonstrates that the fractional rms scatter app(A) in 
the flat field averaged over different size scales A does 
not scale with 1 / ^/A, a clear sign that the flat-fielding er- 
rors are not purely statistical. For scales larger than the 
gaussian FWHM of the seeing disk (~l")j subtle extended 
light from sky objects may not be entirely removed by the 
supersky flat median process, creating an increase in the 
flat-fielding residuals on these scales. On the largest scales, 
systematic errors are nearly an order of magnitude larger 
than those due to counting statistics in the fiat-fields. 

The fractional fiat-fielding uncertainties 5ff(^') from 
Table || must be multiplied by the total unmasked fiux 
in area A' on frame i, and then combined to yield the 
flat-fielding uncertainty within an area A' on the mosaic. 
Since the science frames were dithered by more than 10" in 
each direction, larger than any area considered here, any 
(x, y) position on the mosaic is constructed with images 
that were flat-fielded at different positions on the physical 
detector. Thus the flat-fielding uncertainties of individual 
frames can be treated as being independent and added in 
quadrature. For the mosaic we thus have 



The photon noise is essentially uncorrelated over areas 
larger than the FWHM of the PSF, so that it is given 
by the square root of the number of electrons within that 
area. For a given frame, we thus compute the uncertainty 
due to photon noise within areas Apsp comparable to the 
PSF, and then add these in quadrature. The uncertainties 
for individual frames are independent, and can be added 
in quadrature to yield the total uncertainty due to photon 
noise within an unmasked area A' on the mosaic. Since 
the uncertainties are proportional to the square root of 
the number of electrons but are then added in quadrature, 
for any area A' > Apsp, the resulting photon noise is 



apN{A') = C 



SM 



CsM 



\ 



Nf N^/NpsF NpsF 

E E E 



\ 



EE' 



where N^/Npsp is the number of PSF-sized areas within 
A'. Note that this is nothing more that the square root 
of the total number of electrons recorded by all unmasked 
pixels from all frames contributing to the area A' of the 
mosaic. In all of our work, we will only quote surface pho- 
tometry for A' > Apsp so that all points in our surface 
brightness profiles are independent, and so that the for- 
mulation above can be used to calculate photon noise. 

A. 4. Sky Subtraction: ass 

Sky subtraction introduces the same systematic uncer- 
tainty to every position in the mosaic. The determi- 
nation of the sky values and their uncertainties dS 
(Sr ± SSr = 16651.5±0.4 e- pix-i, 5y ± = 
2950.2 ±0.2 e' 
Sect. 



pix 



were discussed in Sect. 3.4 and 



4.1 



Since the sky values are determined from the 
mosaic itself, the normalization factor Csm is already con- 
tained in these values. We have then simply 



assiA')^N'pSS 



<jpp{A') = C 



SM 



\ 



N, 



N' 



E P^H^oE' 



Csm <7fp{A'] 



\ 



Nf 



N' 



E< 



where j runs over all unmasked pixels in the area A' on 
frame i and aFp{A') can be pulled outside the sum be- 
cause the masking is performed on the mosaic and thus is 
identical for all individual frames i. 



A. 5. Absolute Calibration: ctcal 

The absolute calibration, or transformation of our sur- 
face brightness photometry to a standard system, is not 
of primary importance to many of our scientiflc results 
since relative measurements from one portion of the mo- 
saic to other portions are more relevant. Nevertheless, as 
explained in Sect. 3.5, all absolute measurements have a 
fractional uncertainty of acAL{A') ~ 0.05 due to errors 
in the absolute calibration. Thus, for absolute quantities 
we must also consider 

<Jcal{A') - acAL{A') Fsm{A') 
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Again, because the caUbration is performed on the mosaic 
itself, the normahzation factor Csm is already contained 
in these values. 



A. 6. Mosa icing: au 

If the image of ESO 342-G017 had been formed without 
mosaicing, then the total normalization constant Csm in 
the first equation of the appendix would not be present 
and thus would introduce no uncertainty in the final pho- 
tometry. (The uncertainty in Csm is dominated by the 
uncertainty in ' Si)~^] we ignore here the very much 
smaller uncertainty in ^.) With mosaicing, relative errors 
related to the uncertainty in the quantity (X^f^^ Si)~^ may 
be introduced between different parts of the mosiac. We 
consider, therefore, the mosaicing uncertainty of photom- 
etry in submosaic SM relative to the fiducial submosaic 
SMq (taken to be central submosaic containing flux from 
all ESO 342-G017 frames), to be 



(A')^ v/(5Csm)2 + (5Csmo)2EE' 



C^SM I^SmI'^'^O^ + CsMo SsMo('^'^«)^ V^V^ 

i 3 



c 

where the (SSi)^ are the uncertainties in the median sky 
values of individual frames summed over those frames con- 
tributing to the indicated submosaic. 

A. 7. Surface Briglitness Fluctuations: aL 

In outlining a method for determining extragalactic dis- 
tances, Tonry & Schneider (1988) derive an expression 
for the intrinsic variations in an elliptical galaxy or spi- 
ral galaxy bulge. This fluctuation in surface brightness is 
due to the counting statistics of a finite number of un- 
resolved stars contributing flux to each pixel of a CCD 
image. 

The fluctuations in a single pixel is 



([Tonry fc Schneider 19881 ) 



= 9ti'-^? 10 



-0.4(M-mi) 



where g is the mean number of counts due to the galaxy 
alone (X^T^ 6~ — X^f' '^Oi ^ is the single exposure inte- 
gration time (seconds), d is the distance to the source 
in parsecs, and mi is the magnitude corresponding to 
1 count per pixel per second in the final image. This 
equation, however, does not take into account the effects 
of seeing, which strongly reduce the apparent brightness 
fluctuations. This reduction is significantly greater than 
l/^/n, where n is the total number of pixels, for spatial 



bins larger than the seeing PSF (Morrison et al. 1994). 
Because of the large distance to ESO 342-G017, surface 
brightness fluctuations make a very small contribution to 



our total error budget. Therefore, rather than simulating 
the effects of seeing at various binning scales (as done by 
[Morrison et al. 1994 ), we compute this error as an upper 
limit, and simply scale it to our bin sizes using only a 
l/v^ factor. 

For our data, we take i=600 seconds, fi=102Mpc, 
mi(i?)=21.6mag/sq arcsec and mi(y)=21.3mag/sq arc- 
sec, and M(i?)=-0.7 and M(F)=-0.3 ( [Tonry et al. 19901 ). 



A. 8. The Total Error Budget 

We are now in a position to combine these different sources 
of uncertainty to arrive at an error budget for our sur- 
face brightness photometry of ESO 342-G017. The read 
noise, flat-fielding, photon noise, mosaicing and intrinsic 
surface brightness fluctuation uncertainties are all inde- 
pendent and statistical, and so can be added in quadra- 
ture, so that 



o'Stat(^') — 



'RN 



'FF 



iA') + al.^iA') + aUA') + aliA') 



This is the statistical error that will cause random scat- 
ter in our surface brightness profiles. The sky subtraction 
and absolute calibration uncertainties are systematic in 
that all mcasurmcnts in the mosaic will be affected in the 
same way by errors in these derived quantities. Sky sub- 
traction errors will systematically change the slope of the 
surface brightness profiles; calibration errors will shift the 
profiles by a constant amplitude. All profiles that we dis- 
play in Sect. ^ have statistical error bars computed as de- 
scribed in this appendix. The effect of sky subtraction and 
calibration uncertainties must be considered separately. 
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